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Abstract 

The time-dependent, geometric method for high-energy hmits and inverse 
scattering is applied to nonrelativistic quantum particles in external electro- 
magnetic fields. Both the Schrodinger- and the Pauli equations in and 
are considered. The electrostatic potential Aq shall be short-range, and 
the magnetic field B shall decay faster than |x|~^/^. A natural class of cor- 
responding vector potentials A of medium range is introduced, and the de- 
cay and regularity properties of various gauges are discussed, including the 
transversal gauge, the Coulomb gauge, and the Griesinger vector potentials. 
By a suitable combination of these gauges, B need not be differentiable. The 
scattering operator S is not invariant under the corresponding gauge transfor- 
mations, but experiences an explicit transformation. Both B and Aq are re- 
constructed from an X-ray transform, which is obtained from the high-energy 
limit of S. Here previous results by Arians and Nicoleau are generalized to the 
medium-range situation. In a sequel paper, medium-range vector potentials 
are applied to relativistic scattering. 
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1 Introduction 



The scattering theory of a nonrelativistic quantum particle in an electromagnetic 
field will be discussed under weak decay- and regularity assumptions on the magnetic 
field. Consider first the corresponding classical dynamics, i.e., the Lorentz force: 

mx = e(E + X X B) , (1) 

where m is the mass and e is the charge of the particle, E(x) is electrostatic field 
strength, and B(x) is the magnetostatic field. More precisely, B = /ioH is the 
magnetic fiux density, and H is the magnetic field strength. The field is described 
in terms of a scalar potential Ao(x) and a vector potential A(x) according to E = 
— grad Ao and B = curl A. (Note that there is an alternative system of units of 
measure in use, such that B = H = curl A, where c is the speed of light.) Now 
(1) is equivalent to a Hamiltonian dynamic system with the Hamilton function 

i/(x, p) ■.= T^(p- eA(x))' + eAo(x) , (2) 
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where p = mx + eA(x) is the canonical momentum. A nonrelativistic quantum 
particle is described by a wave function ipi?^) ^ L'^(M.'^, C). Its time evolution is 
determined by the Schrodinger equation ihip = Hip. The self-adjoint Hamiltonian 
H is given by (2), with the canonical momentum operator p = — i^Vx • We set 
h = 1 and e = 1. The Schrodinger operator is describing a spin-0 particle, and the 
similar Pauli operator (53) is describing a particle of spin 1/2, e.g., an electron. If 
Ao{x.) and A(x) decay integrably as |x| oo, i.e., faster than |x|~^, then if is a 
short-range perturbation of Hq := p^/2m. Its time evolution is approximated by 
the free time evolution (generated by Hq) as t ^ ±oo, and the wave operators exist: 

Q±iP:= lim e^^h~^^o^ip. (3) 

f— »itoo 

Here the free state ip is an asymptotic state corresponding to the scattering state 

as t — +00 or t — > — oo, respectively. The scattering operator S := 
is mapping incoming asymptotics to outgoing asymptotics. The vector potential 
A is determined by the magnetic field B only up to a gradient. Under the gauge 
transformation A' = A + grad A, the Hamiltonian H is modified, but the scattering 
operator 5* is invariant if A(x) ^ as |x| ^ oo. 

If A does not decay integrably, the short-range wave operators (3) need not exist. 
The time evolution generated by H may be described asymptotically in terms of 
long-range scattering theory, i.e., by modifying the free time evolution. Loss and 
Thaller [21, 36] have shown that the unmodified wave operators (3) still exist, if 
A(x) is transversal, i.e., x- A(x) = 0, and A(x) = 0(|x|-(i/2+5))_ ^2), (A(x))2 is 
short-range, but A(x) ■ p is formally long-range. It is effectively short-range, since 
A(x) = —X X G(x) with G(x) short-range, and A ■ p = G ■ L with the angular 
momentum L = x x p. This approach generalizes to vector potentials A with the 
following property, which will be called "medium-range" : the transversal component 
of A(x), i.e., orthogonal to x, is 0{\x.\~^^^'^~^^^), and the longitudinal component, i.e., 
parallel to x, is decaying integrably. — The aims of the present paper are: 
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• A class of medium-decay magnetic fields B is considered, such that there is a 
medium-range vector potential A with curl A = B in S', i.e., as a tempered 
distribution. The construction of A requires a decay B(x) = (9(|x|~''^/^+^^) as 
in the case of the transversal gauge, but the local regularity required of B can 
be weakened. The wave operators are obtained analogously to [21, 7]. 

• The decay- and regularity properties of various gauges are discussed, and the 
role of gauge transformations is emphasized: since the scattering operator is 
not invariant in general under the substitution A ^ A' = A + grad A when 
A and A' are medium-range, we must extract gauge-invariant quantities from 
S (only these may be observed in a physical experiment). 

• The corresponding inverse scattering problem can be solved by obtaining the 
X-ray transform of A from the high-energy limit of S. This was done by 
Arians [2] for short-range A under low regularity assumptions, and by Nicoleau 
[24] for smooth A of medium-range. Here these results are extended to low- 
regularity A of medium range. The inverse problem of relativistic scattering 
with medium-range A is addressed in [19], combining the techniques of [22, 36, 
18, 37], including obstacle scattering and the Aharanov-Bohm effect as well. 

Only fields and particles in and M.^ are considered here. A measurable function 
: M'^ ^ M is a scalar potential of short range, if the multiplication operator Ao(x) 
is Kato-small with respect to Hq = ^p"^, and if it satisfies the Enss condition 

/ Ao(x)(i7o + i)"'i^(|x| >r) rfr < cx) , (4) 
Jo 

where F{. . .) denotes multiplication with the characteristic function of the indicated 
region. This condition is satisfied, e.g., if Aq G i^foc ' if decays as |x|~'^ with 
/i > 1. The magnetic field B : M*^ ^ W^' , u' := z/(z/ — l)/2, corresponds to a 2-form. 

Definition 1.1 (Decay Conditions) 

1. Consider a magnetic field B : M'^ ^ or B : ^ M, which is in L^{W) for a 
p> u. It is of medium decay, if it satisfies a decay condition |B(x)| < C|x|~'' for 
a fi > 3/2 and large |x|. In the case ofM."^, we also require that divB = in S' . 

2. A vector potential A : W ^ W is of medium range, if it is continuous and 
satisfies |A(x)| < C|x|~^ for some fi > 1/2. In addition, the longitudinal part 
A(x) • x/|x| shall decay integrably, i.e., 

J sup I |A(x) ■ x|/|x| |x| > r I dr < oo . (5) 

Note that A(x — xq) is of medium range as well. The magnetic fields of medium 
decay form a family of Banach spaces, cf. Cor. 2.5. By the decay at oo, we have 
B G for a p < 2 in addition. The local regularity condition p > u on B enables 
A to be continuous. Since B need not be continuous, the case of an infinitely 
long solenoid is included, where B : ^ R is the characteristic function of the 
solenoid's cross section. The Coulomb gauge vector potential satisfies div A = 0, it 
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is determined uniquely by curl A = B and A(x) ^ as |x| ^ oo. (It is called 
transversal in QED, because p ■ A(p) = 0.) The transversal gauge vector potential 
satisfies x • A(x) = 0, it is determined uniquely by B if it is continuous at x = 0. 
The Griesinger gauge is introduced in Sec. 2.3, motivated by [13]. 

Theorem 1.2 (Vector Potentials) 

1. //A, A' are medium-range vector potentials with curl A' = curl A in S' , then 
there is a C"^ -function A : M'^ ^ M with A' = A + grad A. Moreover, the homoge- 
neous function A(x) := limj.^oo A(rx) exists and is continuous for x 7^ 0. 

2. If B is a magnetic field of medium decay, consider the vector potential A with 
curl A = B m S' , given in the Griesinger gauge. It is of medium range, moreover it 
decays as 0{\x.\~^'^~^'') , ifB decays as (9(|x|~^) and3/2 < ji <2. IfQ is continuous, 
then the transversal gauge vector potential has the same properties. 

3. Suppose that B satisfies a stronger decay condition with fi > 2. Then the Goulomb 
gauge vector potential is of medium range, too. Moreover, A is bounded by C|x|~^ 
in all of these gauges. In M.^, the flux o/B through almost every plane vanishes, and 
the Goulomb vector potential is short-range. In M^, the flux $ o/B is finite, and the 
Goulomb vector potential is short-range, iff'^ = 0. 

4. If B is a magnetic field of medium decay, there is a special choice of a medium- 
range vector potential A = A^ + A*" with curl A = B m S' , where is short-range 
and continuous, and A^ is transversal and , with |c}jA^(x)| < C|x|~'^, /i > 1. 
In addition, divA is continuous and decays integrably. 

Moreover, the Griesinger gauge A and the Coulomb gauge A are regularizing, i.e., all 
diAk have the same local regularity as B. For the transversal gauge A, diA^ exists 
only as a distribution in general if B is continuous. In [36], it is remarked that 
the transversal gauge vector potential is better adapted to scattering theory than 
the Coulomb gauge vector potential, if B(x) = (9(|x|~'^^/^+'^^) and B is sufficiently 
regular. The Coulomb gauge is superior in other cases: 

Remark 1.3 (Advantages of Different Gauges) 

Suppose B : R*^ — s> My' is a magnetic field of medium decay with B(x) = (9(|x|~'^): 

1. If3/2<yU<2, the Coulomb gauge vector potential is not of medium range in 
general, and the wave operators (3) need not exist. The transversal gauge vector 
potential may be used if B is continuous. The Griesinger gauge works in any case. 

2. If/i>2, z/ = 2, and /]g2 B lix 7^ 0, then the Coulomb gauge vector potential is of 
medium range as well, i.e., its longitudinal component is decaying integrably. It is 
preferable to the transversal gauge because of its better local regularity properties, 
and because the Hamiltonian (2) is simplified due to p ■ A = A ■ p. 

3. If > 2 and = 3, or z/ = 2 and /ig2 B da; = 0, the Coulomb gauge vector potential 
is short-range, but the transversal gauge or the Griesinger gauge is short-range only 
in exceptional cases. 

Under the assumptions of item 3, it is natural to use only short-range A, and the 
scattering operator S is gauge- invariant. In the medium-range case, we have a family 
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of scattering operators, which are related by the transformation formula (6): 

Theorem 1.4 (Gauge Transformation, Asymptotics, Inverse Scattering) 

Suppose that B is a magnetic field of medium decay in or M? , A is any medium- 
range vector potential with curl A = B m S' , and Aq is a short-range electrostatic 
potential according to (4). 

1 . For the Schrddinger- or Pauli operator H , the wave operators VL± exist. Consider 
a gauge transformation A' = A + gradA and A(x) = limr_>oo A(rx) according to 
Thm. 1.2, and denote the operators corresponding to A' hy H' , 5". The wave 
operators and scattering operators transform under a change of gauge as 

n'^ = e^^W l]±e-i^(±P) S' = ei^(P) S e'^^^'P^ . (6) 

2. Consider translations in momentum space by u = uuj, G S^~^ . The scattering 
operator S for the corresponding Schrddinger- or Pauli equations in and has 
the asymptotics 

s-lim e~^"^ ^e^"^ = exp {i T oj ■ A(x + u:t) dt} . (7) 

B is reconstructed uniquely from the relative phase of this high- energy limit of S . 
{The absolute phase is not gauge-invariant, thus not observable.) Under stronger 
decay assumptions, error bounds and the reconstruction of Aq are given in Sec. 5.3. 

Item 1 is due to [31] for u = 2. Item 2 will be proved with the time-dependent 
geometric method of Enss and Weder [9]. It is due to Arians [2, 3] for A of short 
range, and in addition for B : — > M of compact support. The inverse scattering 
problem was solved before in [11] for A of exponential decay, and in [24] for C°°-A 
of medium range in the transversal gauge. Analogous results for the Aharanov- 
Bohm effect are discussed by Nicoleau [25] and Weder [37]. The asymptotics for 
C°°-vector potentials of medium range or long range are obtained in [25, 30, 39] 
using the Isozaki-Kitada modification (cf. Sec. 4.5). 

Remark 1.5 (Gauge Invariance) 

1. Gauge freedom has two sides to it: we may choose a convenient gauge to simplify 
a proof, but if the result depends on the gauge, it may be insignificant from a 
physical point of view. Cf. Sec. 6.1. Therefore existence and high-energy limits of 
the wave operators are proved in two steps: first by employing the nice properties of 
the special gauge A = A* + A'' according to item 4 of Thm. 1.2, and then the result 
is transfered to an arbitrary gauge by the transformation (6). Thus (7) is valid in 
any medium-range gauge, and the gauge-invariant relative phase is observable in 
principle. Moreover, this approach shows that only the decay properties of A are 
essential here, while the local regularity properties are for technical convenience. 

2. li u = 2, fi > 2, and the flux of B is not vanishing, it is possible to replace 
the medium-range techniques with short-range techniques plus an adaptive gauge 
transformation, such that the vector potential is decaying integrably in the direction 
of interest. Cf. Cor. 2.8 and [3, 17, 37]. In Sec. 5.3, a different kind of adaptive gauge 
transformation is used, such that the high-energy asymptotics of H are simplified. 
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This paper is organized as follows: Vector potentials are discussed in Sec. 2, including 
the proof of Thm. 1.2. In Sec. 3, B is reconstructed from the X-ray transform of 
A. The direct problem of nonrelativistic scattering theory is addressed in Sec. 4. 
Existence of the wave operators is proved in detail, because the same techniques 
are needed later for the high-energy limit, but the reader is referred to [21, 7, 1] for 
asymptotic completeness. Sec. 5 is dedicated to the inverse problem, and concluding 
remarks on gauge invar iance and on inverse scattering are given in Sec. 6. 
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2 Fields and Gauges 

To construct medium-range vector potentials, we are employing vector analysis 
on and under low regularity assumptions, controlling the decay at infin- 
ity. Some references to similar results for Sobolev spaces over domains in are 
included, and [29] is a standard reference for vector analysis on manifolds using 
distributional derivatives. The following notation will be employed in the case of 
M^. It is motivated by identifying vectors and scalars in with vectors in M^, 
\ = {vi, ^2)*'^ ^ {vi , V2 , 0)*'' and w ^ (0 , , w)*'': 

X X V := X1V2 — X2V1 curl V := diV2 — d2Vi 

X X w := {x2W , —XiwY'^ curlw := {8210 , —diwY'^ . 



2.1 Gauge Transformation of A 

Suppose A : M'' ^ M'' is a medium-range vector potential according to Def. 1.1. 
Assume curl A = in S', i.e., / A-curl rfx = for G ^(M^ M) or (/> e 5(M^ M^), 
respectively. To show that A is a gradient, one can mollify A and employ a density 
argument, or apply mollifiers to the line integral over a closed curve. Here we shall 
use test functions, and define A : M*^ — ^ M by the Poincare formula for closed 1-forms, 

A(x) := f\-A{sx.)ds . (8) 
Jo 

For (j) G iSflR'^, M) and a unit vector a; G S'^~^ consider 

A(x)lc;- V0(x)cia; (9) 
/ / X- A(sx)cj- V0(x)c?srfx (10) 

JE" JO 







A(sx) ■ curl (o; X x</)(x)) dsrfx — (11) 
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-/ /^a;-A(sx)((z/-l)0(x)+x- V0(x))dsrfx. (12) 

JR" Jo ^ ' 

This identity is verified with a x (b x c) = (a ■ c)b — (a • b)c. Now (11) is vanishing 
by curl A = in S' . The substitution (s, x) (s, y) with y = sx in (12) gives 

- / a; . A(y) (iv - \)s~^ <\>{y j s) + s"^- V ■ Vx0(y/s)) ds dy (13) 

JM" Jo ^ ' 

a;.A(y)p->(y/.)l\ = / a; • A(y) 0(y) d?/ . (14) 



This shows A = grad A in 5', and we have A G since A is continuous, thus 
A = grad A pointwise. Moreover, hne integrals of A are path- independent. Now 
consider 

A(x) := lim A(rx) = lim / x • A(tx) dt . (15) 

Since ^ ■ A(x) is short-range, convergence is uniform for |x| > i?, thus A is con- 
tinuous on W \ {0}. (A is 0- homogeneous, and discontinuous at x = unless it is 
constant.) If A is short-range, it is well-known that A can be redefined such that 
lim|x|^oo A(x) = 0. This follows from the estimate A(x) = o(|x|~^), see, e.g., [17, 
Lemma 2.12]. Item 1 of Thm. 1.2 is proved. ■ 

Remark 2.1 (Poincare Lemma I) 

1. The same proof shows the following version of the Poincare Lemma: Suppose 
1] C M^ p > z/ > 2, and A e L\^X^, W) with ^ A^Sfc^ - A^di4)dx = for all 
G C^{Q). If Q is starlike around x = 0, define A : ^ M by (8) a.e., then A is 
weakly differentiable with grad A = A almost everywhere. 

2. If Q is simply connected, then A is obtained piecewise. For general Q, A need 
not exist globally. Now A is a gradient, iff /q A • vrfx = for all v G C^{^1, W) 
with divv = 0. Under this assumption on A G L^^^^Q, M.'^), A is constructed by 
employing mollifiers [32, 12]. For suitable Q, this result implies the Helmholtz-Weyl 
decomposition of L^(f2, R'^), p > 1, which is important, e.g., in fluid mechanics. 



2.2 The Transversal Gauge 

Suppose that B : ^ B : _^ M is a magnetic field of medium decay 

according to Def. 1.1. The transversal gauge vector potential A : M'^ — >■ is defined 
a.e. by the Poincare formula for closed 2-forms, 

A(x) := -X X /\b(sx) rfs . (16) 
Jo 

Proposition 2.2 (Transversal Gauge) 

Suppose that B : M'^ ^ M*^' is a magnetic field of medium decay with B(x) = 
(!?(|x|~^) as |x| oo, fi > 3/2. In the transversal gauge, the vector potential A is 
defined by (16). Assume in addition that B is continuous. Then 

1. A is continuous and satisfies A(x) = (Ddxl^-*^) if fi> 2, A(x) = (!?(|x|~^ log |x|) 
if H = 2, and A(x) = 0{\:>c\~^''~^^) if fi< 2. Since x ■ A(x) = and fi > 3/2, A is 
of medium range. 
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2. We have curl A = B m S' , but the weak partial derivatives diAk and div A do 
not exist in general. 



B is required to be continuous, to ensure that A is continuous. (More generally, 
B may have local singularities c|x — xo|~^^~'''', or a jump discontinuity on a strictly 
convex line or surface, but a jump discontinuity on a line through the origin is not 
permitted.) The decay properties are given in [21, 22, 36, 7] for 3/2 < /i < 2, and 
in [25] for /i 7^ 2. To achieve that A G with all derivatives decaying integrably, 
we would have to assume B G C"*^ with derivatives decaying faster than |x|~^. 

Proof: 1. Define b{r) := sup|x|=r |B(x)| for r > 0, then b{r) = 0{r~^) as r — 
Now |A(x)| < \^\-^ fl,'"^ r b{r) dr gives the desired estimates (which are sharp). 
2. Consider z/ = 3 and a test function (p G 5(]R^, M^). We have 



00. 



+ 



Afx) ■ curl0(x) dx 

X X sB(sx)) ■ curl 0(x) ds dx 
/)(x) I ds dx 



f /\sB(sx)-Vfx 
-f rsB(sx)- (0(x) + (x-V)0(x))rfs(ix 



(17) 
(18) 
(19) 
(20) 



since (a x b) ■ (c x d) = (a ■ c)(b ■ d) — (a ■ d)(b ■ c). The integral (19) is vanishing 
because of divB = in S'. In (20) we substitute y = sx and obtain 



B(y) ■ (s^'^{y/s) + s-%y ■ Vx)0(y/s) dsdy 



s=0+ 



dy 



B(y)-0(y) dy. 



B(y) ■ [s~'<Piy/s 
In dimension u = 2, we have G iS(M^, M), and (18) equals 
— / / sB(sx) X ■ V0(x) (is cia; 

™'B(y)s-V-Vx0(y/s)rf5d|/ = / B(y)U(y/s"' 



(21) 
(22) 

(23) 
dy . (24) 

Thus curl A = B in iS'(]R^, M^) or iS'(]R^, M), respectively. — For ly = 2, suppose 
that B(r cos6', rsin6') = (l+r)~^/(6'), where / is singular continuous. Then none of 
the weak derivatives diAk or div A exists in Lj'Qj,(M^). In the following example, the 
derivatives exist but they are not short-range: B(rcos6', rsin^) = (l+r)"^ cos(r^6'). 
Similar examples are constructed in M^. (The condition divB = is satisfied, e.g., 
by B(x) = X X grad5f(x).) H 

Remark 2.3 (Poincare Lemma II) 

1. Suppose ncR^,p> 3/2, and B G Lf^^i^, R^) with ^B ■ grad0c/x = for all 
G C^{Q). If Q is starlike around x = 0, define A(x) by (16) a.e., then the same 
proof shows a ■ curl (pdx = J^B ■ (pdx for G C^{Q, M.^). This vector potential 
is not weakly different iable in general. 
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2. On an arbitrary domain Q, a vector potential A exists if ■ (pdx = for 
all (p e C^{^, M^) with curl0 = 0. If can be replaced with C°°, then A can 
be chosen such that it vanishes at the (regular) boundary dQ: A is constructed in 
[38] by potential theory, and in [13] by a mollified version of (16), see below. These 
vector potentials are weakly differentiable with ||9tv4fc||p < Cp||B||p, 1 < p < oo. 

2.3 The Griesinger Gauge 

For a magnetic field B : M'' ^ M'^' of medium decay, the Griesinger gauge vector 
potential shall be defined by employing a mollifier h G C^(M'^, M) with J^^ hdx = 1: 

A(x) := - f f\{z) - z) X sB{sx + {1 - s)z) ds dz (25) 

JR"^ Jo 

= - - - y)) t'~\t - 1) - y) X B{y) dt dy . (26) 

JR" Jl 

(25) looks like a mollified version of the transversal gauge (16), which is recovered 
formally for h{z) — > 6{z). Note that B is averaged over a ball of radius ^ (1 — s), 
which is shrinking to a point as s — > 1 in (25), or y — > x in (26). It turns out that 
the integral kernel in the latter equation is weakly singular. 

This definition is the exterior domain analog to the construction found by Griesinger 
for interior domains [13], i.e., ds was replaced with — ds. (Her original con- 
struction is less suitable for magnetic fields of medium decay, because it would 
require B(x) = Odxl^^""^), and vanishing flux in the case of M?.) The technique 
goes back to Bogovskij's solution of div v = / [4, 5, 12]. In the case of a bounded do- 
main, the vector field satisfies < Cp||/||p or ||5iAfc||p < Cp||B||p, respectively 
(1 < p < oo). — Item 2 of Thm. 1.2 is contained in the following 

Proposition 2.4 (Griesinger Gauge) 

Fix h G C^{M.'^, M) with /jj^ hdx = 1. Suppose B is a magnetic field of medium 
decay with B(x) = (9(|x|~^), > 3/2. The Griesinger gauge vector potential A is 
defined by (25). 

1. A is continuous and satisfies A(x) = (9(|x|~^) if fi > 2, A(x) = (9(|x|~^ log |x|) 
if H = 2, and A(x) = 0{\x.\~^^~^^) if n < 2. The longitudinal component of A is 
short-range, thus A is of medium range. 

2. A has weak partial derivatives in Li^^{My), and curl A = B almost everywhere. 
But the weak derivatives do not decay as 0{\x.\~^) in general. 

Proof: 1. Choose p > u with B G Lp{W, W'), q := 1/(1 - 1/p), and fix i? > 
such that h{z) = for |z| > R and |B(x)| < c(l + |x|)^^ for |x| > R. By Holder's 
inequality we have 

|A(x)| < (|x| + R) \\h\U t s||B(sx + (1 - s)z) x(|z| < i?)||p ds . (27) 

JO 

Note that the norm of B is considered on the ball of radius (1 — s)R around sx. For 
|x| < 2R consider || . . . ||p < ||z ^ B(sx+ (1 - s)z)||p = (1 - s)~'^/P||B||p , thus 

|A(x)| < 3i?||/i|L||B|L (\s{l - 3)-"'^ ds (28) 

Jo 
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is bounded. For |x| > 2/?, the s-interval in (27) is split: 
IR 

a) For < s < — we have 1 - s > 1/3, thus || . . . ||p < S^^/fUBHp =: ci . 

|x| + R 

IR 

b) For ■; — ; < s < 1 we have |sx+ (1 — s)z| > R and |sx+ (1 — s)z| > s|x|/2, 

|x| + ~ ~ ' ^ ^ ' ~ ' ^ y I _ I 1/ , 

therefore || . • • ||p < C2(l + s|x|)~^. 

Now ||B(sx+ (1 - s)z)x(|z| < R)\\p < C3(l + s|x|)-'' for < s < 1, and (27) 
yields the desired estimate for |A(x)|, which is optimal. The stronger bound for the 
longitudinal component is obtained by replacing the factor (|x| + i?) ^ i? in (27). 

Write (26) as A(x) = /jg^ ^(x, x — y) B(y) dy. The kernel (^(x, z) is bounded by 
c(l + Ixl*^^^) |z|~^'^~-^\ analogously to (38) below. Thus it is weakly singular, and 
B G implies that A is continuous (adapting Thm. II. 9. 2 in [12]). 

2. We show curl A = B in 5' by the techniques from Sec. 2.2, and establish the 
higher regularity of the Griesinger gauge afterwards (which was not possible for the 
transversal gauge). Consider z/ = 3 and cj) G iS(M^, M?). As in (19) and (20) we have 



+ / A(x) ■ curl0(x)dx (29) 
= 11 f\(z)sB(sx+(l-s)z)-VJ(x-z)-(f)(x))dsdzdx (30) 

- [ [ /\(z)sB(sx+(l-s)z)- f0(x) + ((x-z)-Vx)0(x))c/srf2rfx .(31) 

In (30), the x-integral is vanishing for a.e. s and z, because divB = in S'. In (31), 
the substitution (s, z, x) (s, z, y) with y = sx + (1 — s)z yields 

^\(z)B(y)- 



s-V(z + (y - z)/s) + s~%{y - z) ■ Vx)0(z + (y - z)/s)) ds dzdy (32) 
/ h{z) B(y) ■ [s- V(z + (y - z)/s)l ' dz dy (33) 
h{z)dz / B(y)-0(y)rf2/, (34) 

thus curl A = B in 5'(]R^ M^). For z/ = 2 we obtain [(j){z + (y - z)/s)]^_^ = (f){y) 
analogously to (24). (The same technique works for ds, but for u = 2 we have 
[0(z + (y - z)/s)]]^^_ = 0(y) - 0(z), thus curl A(x) = B(x) - /i(x) J^, B(y) rfy.) 

The existence of all weak derivatives diA^, 1 < i, k < u, is shown first for B G 
(without restriction on divB), and then a density argument covers the general case. 
To simplify the notation, we consider only div A: (26) implies, as a principal value, 

divA(x) = / K(x, x-y) ■B(y)rfy with (35) 

Vh{^-tz)t''~\t-lfdt (36) 
/ VhU-ri^\r''~\r -\z\fdr . (37) 

J\z\ ^ \z\' 
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The most singular contribution is 

7 roo / 7 \ / \ 

~^P+i^i) '^^(^-^^)^'^^ = ^((i + i^r)|2r') • (38) 

The Calderon-Zygmund Theorem [6], cf. [13, 12], shows that the principal value 
integral (35) is well-defined, and 



divAllp < Cp||B||p . (39) 



Approximating the given B G with B G C^, (39) shows div A G L^^^ . 

Now diAk G L^Q^iW^) is shown analogously, and curl A = B in iS' implies curl A = 
B a.e.. 1 do not know if diAk G Lp(M'^), or if these weak derivatives decay integrably 
as |x| — > oo (short-range terms). But assuming that all (9j/lfc(x) = C(|x|~'^), i.e., 
having the same decay as B(x), would imply A(x) = 0{\x.\~^^~^^), which is not true 
in general if yU > 2. ■ 

Corollary 2.5 (Banach Spaces) 

For R > 0, p > u, and fi > 3/2, define a norm 

\\B\\n,,,, := ||B(x)x(|x| < R)\\p+\\ |xrB(x)x(|x| > i?)|U , (40) 

and denote hy M-R^p^^ the Banach space of magnetic fields with finite norm, and with 
divB = in S' if u = 3. The vector space of magnetic fields with medium decay 
is the union of these Banach spaces. For fixed h G , the proof of item 1 above 
shows that the Griesinger gauge is a hounded operator Jvlji^^^ — > C°(M'^, W). 



2.4 The Coulomb Gauge 

In the Coulomb gauge, the vector potential A is defined by 

A(x) := / X B{y)dy (41) 

uju JR" |x — y| 

with := I'S'-^I = 27r and uj-^ := l^^l = 47r. 

Proposition 2.6 (Coulomb Gauge) 

Suppose B : M'^ ^ W^' is a magnetic field of medium decay with B(x) = (9(|x|~'^), 
/i > 3/2. The Coulomb gauge vector potential A -.W is defined by (41). 

1. A is continuous and weakly differentiable with \\diAk\\2 < ||B||2 . It satisfies 
div A = and curl A = B a.e.. 

2. If ^> 2, then A is of medium range. For v = 3 it is short-range. For v = 2 it is 
short-range, iff the flux is vanishing: /]g2 B cix = 0. 

3. 7/3/2 < /i < 2, then A satisfies A(x) = 0(|x|-(^-i)), or A(x) = C(|x|-i log |x|) 
if fi = u = 2. But the longitudinal component of A does not decay integrably in 
general, and A need not be of medium range. 
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Thus the Coulomb A has better differentiabihty properties than the transversal 
gauge, even if B is continuous. The decay properties are better for /i > 2 (and 
vanishing flux), but not sufficient in general if yU < 2. — In [23], |A(x)| is estimated 
in and in exterior domains. We shall employ a convolution of Riesz potentials: 

Lemma 2.7 (Riesz Potentials) 

For z/ G N and < a, 13 < u with a + (3 > v , we have the convolution on 

|x|-" * |x|-^ = C„,^;^|x|-("+^-^) . (42) 

This formula is obtained form an elementary scaling argument and convergence 
proof. The constant is determined in [14, p. 136]. 

Proof of Prop. 2.6: 1. A is continuous since the integral kernel in (41) is weakly 
singular and B G for some p > u, cf. Thm. II. 9. 2 in [12] (the condition of a 
bounded domain is overcome by splitting B). The convolution (41) is obtained 
by differentiating the fundamental solution of the Laplacian: we have A = curl U 
with — AU = B. This implies divA = and curl A = B in S'. The Fourier 
transforms satisfy ip ■ A(p) = and ip x A(p) = B(p), thus p^A(p) = ip x B(p). 
It remains to show that A is weakly differentiable. Now |9jy4fc(p)| < |B(p)| a.e., 
thus ||(9jAfc||2 < ||B||2. We have B G for pi < p < P2 , with pi < 2 and p2 > v. 
By the Calderon-Zygmund Theorem [6, 13, 12], ||9jAi.||p < Cp||B||p for pi < p < P2 ■ 

2. Now /i > 2, and we may assume 2 < yU < 3. Split B = B*^-"^) + B^^^ , such that 
(z jj> hag compact support, with B*^^) = in the case of M^, and such that 
|B(2)(x)| < c|x|-^ for X G (divB^ = is not required). Split A = A^^) + A^^) 
according to (41), i.e., by applying the convolution to B'-*'' individually. If z/ = 3, 
we have |A'^^)(x)| = C(|x|~^), since |x — y|~^ = C(|x|~^) for y G supp(B'^-^)), and 
|A*^^''(x)| = (9(|x|~'^^~^'') by Lemma 2.7, thus A is short-range. For z/ = 2 we claim 



A(x) = ' \ ( B(y) dy + 0(|x|-(^-i)) . (43) 




The integral kernel of (41) is decomposed as follows: 
x-y X (xx y) X (x-y) - (x-y)(x-y) 



O 



I |2 Ii2' Ii2i |2 Ii2 

|x — y| |x| |x| |x — y| |x| ^ |x| |x 

When applying this kernel to B*^^-* , the first integral is vanishing and the second is 
(9(|x|~^). Applying it to B*^^-* , the first integral gives the leading term in (43), and 
the second integral is bounded by 

^dy=^f |x-y|-i|y|-(^-^'rfy = c'|x|~('^-^) (44) 



X X — V X 



by Lemma 2.7. Thus (43) is proved. (For B of compact support this is due to [35].) 
If the flux of B is vanishing, then A is short-range. If not, then A is still of medium 
range, since the leading term is transversal. 

3. Now 3/2 < /i < 2. Split B and A as in the previous item, then A^^^(x) = 
0(|x|^2)^ |B(2)(x)| < clxl^'', and Lemma 2.7 gives |A(2)(x)| < c'|x|~(^~^), except 
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for the case of /i = = 2: then compute the bound exphcitly for c(l + |xp) 
Consider the vector potential A(x) given by 



1 / xi{xl - xl + 1) 



or 



f xi{xl - xl + xl + 1) \ 

X2[x\ — x\ — x\ — \) 

\ / 



respectively, and define B := curl A. We have B G with B(x) = 0(|x| ^) and 
divA = 0, thus A is the Coulomb gauge vector potential for the medium-decay 

X\ X\ 

B. Now 0) = —2 or 0, 0) = , respectively, shows that the 

X]^ + 1 X;^ + 1 

longitudinal component is not short-range, and A is not of medium range. ■ 

In the last example, A(x) := limj.^oo A(rx) according to item 1 of Thm. 1.2 does 
not exist for the gauge transformation from transversal gauge to Coulomb gauge. 
The scattering operator exists for the transversal gauge, but not for the Coulomb 
gauge. Note that for rotationally symmetric B : M, the Coulomb gauge agrees 

with the transversal gauge. Therefore, this vector potential combines the regularity 
properties of the Coulomb gauge with the decay properties of the transversal gauge: 
it is weakly different iable and of medium range for all /i > 3/2. 

If /i > 2, the flux of B is finite if z/ = 2, and for = 3, the flux through almost 
every plane is vanishing. When v = 2 and l^^'Qdx ^ 0, consider the natural class 
of medium-range vector potentials satisfying A(x) = (9(|x|~^). The corresponding 
gauge transformations A(x) have the property that A(x) is Lipschitz continuous for 
|x| > e. This class contains the Coulomb gauge, and the transversal gauge as well 
if B is continuous. 

Corollary 2.8 (Adaptive Gauges) 

Suppose V = 2 and B is of medium decay with B(x) = (9(|x|~^) for a fi > 2, 
and $ := /jg2B(ix ^ 0. For any direction E there is a vector potential 
A'^ of medium range, such that A'^(x) = Odxl"-"^), and A'^(x) decays integrably 
as 0{\x.\~^^~^^) within sectors around ±uj. Moreover, divA'^ is continuous with 
divA'^(x) = C(|x|-2). 

Proof: Denote the Coulomb gauge vector potential by A and observe (43). Choose 
a 27r-periodic function / G C^(M, M), such that f'{9) = ^/{2tt) for 9 in intervals 
around arg(±a;). Choose A G C^(M^, M) with A(x) = /(argx) for large |x|. Then 
consider A'^ := A — grad A. ■ 

Similar constructions are found, e.g., in [37]. If B has compact support, A'^ can 
be chosen to vanish in these sectors for large |x|: this is achieved by subtracting a 
gradient, or by a shifted transversal gauge if B is continuous [3, 17, 37]. 



2.5 The Hormander Decomposition of A 

The following lemma is a special case of [15, Lemma 3.3], which is a standard tool 
in long-range scattering theory. (It is used to improve the decay of derivatives of 
the long-range part Aq , where Aq = Aq + Aq .) 
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Lemma 2.9 (Hormander Decomposition) 

Suppose that V E C\W) and V(x) = O(|x|-"^o), mv{x) = C(|x|-'"i) for \-f\ = 1, 
with mo > mi — 1 > 0. For < A < min(l, rrii — 1) there is a decomposition 
V = Vi + V2 such that: Vi E is a short-range potential with Vl(x) = C(|x|~'^) for 
A = max(mo, mi — A) > 1, and V2 G C°° satisfies d'^V2{'x) = Cdxp^j) for I7I = j , 
J G No . Here = rriQ and m'- = max(mo + jA, mi + (j — 1) A) for j G N. 

For a given medium-decay B we want to obtain a corresponding medium-range 
A = A'^ + A'', sucli tliat A* is short-range, A'' is transversal with short-range 
derivatives, and divA is short-range. The transversal gauge does not satisfy our 
requirements because it need not be differentiable, and the Coulomb gauge is weakly 
differentiable and satisfies the condition on div A, but its longitudinal part need not 
decay integrably. For the Griesinger gauge, I do not know how to control the decay 
of the derivatives. Lemma 2.9 cannot be applied directly to B or to a known A, 
because it requires B G or A G C^, respectively. 

Proof of Thm. 1.2, item 4: Suppose B(x) = 0(|x|~'*) for |x| 00 and consider 
the Coulomb gauge vector potential A'^ according to (41). If yU > 2, we may take 
A* := A^ and A'^' := 0, except in the case of = 2 and j"Bdx ^ 0: then A*^' 
equals the first term in (43) for large |x|. Thus we may assume 3/2 < /i < 2, or 
fi = 3/2 + 35 with < 6 < 1/6. A*^ is continuous, weakly differentiable, and satisfies 
A=(x) = C(|x|-i/2-35)^ Choose a function r] G C~(R'^, M) with /^(x) = for |x| < 1 
and ?7(x) = 1 for |x| > 2, and consider the decomposition A"^ = A^^^ + A°° : 

A(°Hx) = -- / {l-vi^-y))r^^ xBiy)dy, 
ujy Jw |x — y| 

A~(x) = -- / r/(x-y)/^i^ X B(y)rfy. 
ujy Jw |x — y| 

Now A(°) is short-range and A°° G C°° with A~(x) = C(|x|~^/2-35)^ r^^ie deriva- 
tives are given by convolutions diA^ = J^Kiki * Bi, where Kiki(x.) = Cdxl"*^) as 
|x| 00 and i^'»fci(x) = for |x| < 1. Thus |i^ifc«(x)| < clxj-^'^-^) for x G W. Split 
B = B(^) + B(2), such that B^^^ has compact support and |B(2)(x)| < c'|x|-2/2-35_ 
The estimate diA'^{x.) = (9(|x|~^/^~^^) is obtained from Lemma 2.7. A medium- 
range A can be constructed as A*-*^^ plus the transversal gauge for curl A°°, but the 
derivatives of the latter need not decay integrably. 

Lemma 2.9 with mo = 1/2+35, mi = 3/2+25 and A = 1/2+5 yields a decomposition 
A°° = A(^) + A(2), such that A^^) is short-range and ^^^^^(x) = Cdxj-'^i) for 
1^1 = j, with m[ = mi = 3/2 + 25 and m'^ = 2 + 35. Now define A'' := A(°) + 
A^^^ = A^ — A(^\ then A^ is continuous, weakly differentiable, and short-range. The 
longitudinal part of A*^^) need not decay integrably. Consider the decomposition 
B = B** + B'' with B^ = curl A"* in S' and B"" = curl A^^^ in C°°, which satisfies 
divB^ = divB^' = 0. We have B^'(x) = C»(|x|-3/2-2'5) ^^^j ^^^^(x) = 0(|x|-2-35). 
Define A*" as the transversal gauge vector potential belonging to B*^', then A*" G 
with A''(x) = (9(|x|^^/^^^^). By differentiating (16) under the integral and an 
analogous estimate, diAl(x.) = Odxl^^"^"^) is obtained. Now A := A* + A*" yields 
the desired gauge. Note that div A = div A^ — div A^^^^ G C°° decays integrably. ■ 
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3 Inversion of X-Ray Transforms 



From the asymptotics of the scattering operator S, we will know the line integral of 
A along all straight lines in M.'^, up to adding a function of the direction uj G S"'^: 

/oo 
UJ ■ A{x + (jjt) dt . (45) 
-oo 

Proposition 3.1 (Inversion of X-Ray Transforms) 

Suppose A is an unknown vector potential of medium range, and the line integral 
(45) is given for all x G M'^ and to G 5''^""'^, up to adding a function f{uj)- Then 

1. The distribution B := curl A G S' is determined uniquely. 

2. Assume in addition that B is a magnetic field of medium decay. On a.e. plane, 
the X-ray transform of the normal component of B is obtained from (52) below. 

Under stronger decay assumptions on A, item 1 is due to [17, 18], and under stronger 
regularity assumptions, item 2 is found, e.g., in [16, 25]. 

Proof: 1. Given G S{W, W'), we will need a vector field iJj G C°°(M^, W) with 

curl0(^) = 2/ yy^(^-y) ^^. (46) 
JR" |y| 

By Fourier transformation, this equation is equivalent to 

ip X 0(p) = ^ (|Pr^(p) - P P ■ ^(P)) (47) 

for some c^, > 0. Choose the following solution G C°°: 

^(p) := ic;^|p|p X 0(p) = c;^[p[~^ ([p[=^ip X 0(p 
^(x) = -c'^ [ [x- u["('^~^) Acurl0(u)rfM . (4. 



To determine the decay properties of ip, split [x[~'^'^~^) = /i(x) + /2(x) such that 
/i has compact support and /2 G C°°, and split = i^i + according to (48). 
Integrating by parts three times shows that ip2 is the convolution of with a C°°- 
kernel, which is 0{\x.\~^'^'^'^^), thus in L^. Therefore ip G L^. Now suppose that 
a{uj, x) = UJ ■ A(x + ujt) dt is known up to a constant depending on uj, then 

/ a{uj, X.) UJ ■ ip{x.) dx = / / UJ ■ A{-K + ujt) UJ ■ ip{-K) dxdt (49) 

JR"^ J-oo 

is known: the unknown constant is canceled, since ip{0) = gives f^uipdx = 0. 
Consider polar coordinates y = ujt, dy = \'y\'^~^du} dt to obtain 

/ / a{uj, X.) UJ ■ ip{-x) dxdu = 2 / / y A(x + y) y ip{'x) ^^^^ 
Js"-^ JR" Jr" jr" ly[ 



JR" J my [y[ ^ jr" 

by (46), thus the distribution B = curl A G S' has been computed. 
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2. Fix unit vectors u;, a; G 5^ ^. If A G C^, thus B is continuous, we have 

/CO POO 
ui ■ + oju + ujt) dt = / {uj x u) ■ 'B((K + U3U + (jJt) dt . (50) 
-oo J —oo 



d 
du 



(Proof by Stokes' Theorem, or with (cDxa;)-(Vx A) = (a;- V)(ci;- A) — (aj-V)(cD-A).) 
Thus the X-ray transform of the component of B in the direction of u3 x a; is obtained 
on every plane normal to a) x cj. (It is natural, but not required, to assume uj-u) = 0.) 
The left hand side of (50) does not depend on the gauge of A, since duf{uj) = 0. 
Now B G L^, and A is of medium range. Then (50) remains true for a.e. x G M*^ 
and a.e. n G M. The proof is given for z/ = 3, a; = (0, 1, 0)*'' and oj = (1, 0, O)*'": 
We have 

/oo roo ru roo 

A2{u, t, xs) dt - A2{0, t, X3) dt = / Bsiv, t, X3)dtdv (51) 
-00 J —00 Jo J —00 

for every u and almost every X3 , by integrating with respect to X3 , an approximation 
argument, and Fubini. For almost every X3 , both sides are well-defined, and the 
right hand side is weakly differentiable with respect to u. (50) may be rewritten as 

/oo 
(cD X a;) ■ B(x + a;t)rft , (52) 
-00 

where us ■ V denotes a weak directional derivative, which exists for a.e. x. ■ 

B is reconstructed from the X-ray transform according to [14, 9, 18]. We have 
considered the normal component of B on almost every plane: 

Remark 3.2 (Trace of B on a Plane) 

Suppose B : ^ is a magnetic field of medium decay. If B is continuous 
except for a jump discontinuity transversal to a surface, the condition divB = in 
S' implies that the normal component of B is continuous. If B G with p > 3, I do 
not know if there is a kind of trace operator, which defines the restriction of a; • B to 
every plane normal to a;, in Lj'Q^(]R^) or in L'^iM?). (This restriction is well-defined 
as a distribution in iS'(]R^, M), by employing a vector potential.) 



4 The Direct Problem of Scattering Theory 

Now vector potentials of medium range are applied to a nonrelativistic scattering 
problem: this section contains the proof of Thm. 1.4, item 1. 



4.1 Definition of Hamiltonians 

Our Hilbert space is 7i = L'^iW, C) for the Schrodinger equation, and 7i = 
L'^{W , C^) for the Pauli equation. In the latter case, the Pauli matrices CTj G C^^^ 
are employed [36]. The free time evolution is generated by the free Hamiltonian 
-^0 ~ ~2m^ ~ 2mP^- self-adjoiut with domain = W^, a Sobolev space. In 



16 



an external electromagnetic field, the Pauli Hamiltonian is defined formally by the 
following expressions: 



H = -^kp-A)2-o-.Bl+Ao (53) 
2m L -I 

= — [p^ - 2A ■ p + i div A + - cr ■ bI + Ao . (54) 
2m L J 

Domains will be specified below by employing perturbation theory of operators and 
quadratic forms, cf. [28, 27]. The following quadratic form is needed as well: 



1 



2m 



2 r . , 1 



qAi^,^):= — {p-A)i> + U, (_ o--B + Ao)^ . (55) 



2m 



Lemma 4.1 (Pauli- and Schrodinger Operators) 

Suppose B is a magnetic field of medium decay and A is a corresponding vector 
potential of medium range, and Aq is a scalar potential of short range, cf. (4) . The 
Pauli operator H is defined in item 1 : 

1. There is a unique self-adjoint operator H , such that its form domain is the Sobolev 
space , and the quadratic form corresponding to H equals q^ according to (55). 

2. Suppose A' = A + grad A is of medium range as well, and define H' in terms of 
qA' analogously to item 1. Then Dh> = e^^^^^Z)// and H' = e^^(^)i/e~^^(^) . 

3. Suppose the distribution divA E S' is a bounded function. Then Dh = W^, and 
H satisfies (54). 

For the Schrodinger operator, the term —a ■ B is omitted, and curl A E S' need not 
be a function. 

Proof: 1. The quadratic form (55) is well-defined on , which is the form domain 
of Hq . It satisfies |5a('0) V')! ^ (^(V^) H^ip) + 0{\\iIj\\'^) for some a < 1: the bound for 
is < 1 since Aq is short-range, and the bounds of the other terms are arbitrarily 
small. Now H is obtained from the KLMN Theorem. 

2. grad A = A' — A is bounded and continuous. The mapping ijj i— > ^"^ijj is unitary 
in Ti and sending to itself. On we have 

gi^W _ A(x)) e^i^ W = p - grad A(x) - A(x) = p - A'(x) , (56) 

thus qA'ii^, = '?A(e~^'^'?A, e~^'^ip) for i/j G W-^. Since H and H' are determined 
uniquely by the quadratic forms, we have H' = e^'^He~^'^ . 

3. Now a sum of distributions is a bounded function. Denoting it by divA, we 
have J^v A ■ grade/) -|- [div A)(f)dx = for (j) E S. (If the distributions curl A and 
div A are in Lf^^ , then all diA^ G Lj^^, by elliptic regularity. But here we include the 
more general case of a Schrodinger operator as well, where div A is a function but 
curl A is not.) The expression (54) is a symmetric operator on S. By the Kato- 
Rellich Theorem, its closure is a self-adjoint operator H with Dj^ = Dhq = W^. We 
compute qA{'4'i = ("^^ Hip) for ip E S. By the KLMN Theorem, iS is a form core 
for H, thus H = H. u 



17 



4.2 Existence of Wave Operators 



Consider a Pauli operator H according to Lemma 4.1. We assume that A is given 
by item 4 of Thm. 1.2, i.e., A = A'' + A^ with A'^ short-range, A^ transversal, and 
divA short-range. Thus H is given by (54). Existence of the wave operators will 
be shown by estimating the Cook integral 

nitoo 

n±ip:= lim e^^^e~^^or^ = ^ + i 



ei^^ - Ho) e-'^ot ^ _ (57) 



The integral is well-defined on a finite t-interval ii & D = Dh = W^. We have 
H-Hq = -—A''-p + V\ V^^ := — |-2A^-p + idivA + A2-<T-B| + Ao . (58) 

777, 9rn ^. J 



Fix V e 5 with i) e C^{W \ {0}) and e > with V^(p) = for |p| < em. Choose 



with 5'(p) = 1 on supp(?/') and consider the decomposition 



< 
+ 



Aogip)Fi\^\>e\t\) 



Aogip) ■ F(|x| <£|t|)e-i^o^V' 



(59) 
(60) 



The first summand has an integrable bound since Aq is a short-range potential (in 
(4), the resolvent may be replaced with g{p))- The second term can be estimated by 
any inverse power of |t| , by a standard nonstationary phase estimate for propagation 
into the classically forbidden region [36], since the speed is bounded below by e. The 
remaining terms in are treated analogously, observing the decay properties of A*, 
A^, div A, and B. In the term A* • p the space decomposition is introduced between 
A* and p. The term A'' ■ p in (58) is controlled with the technique of [21, 36]: write 
A''(x) = —X X G(x) and note that (G x x) ■ p = G ■ L with the angular momentum 
L = X X p. Now L is commuting with Hq , and G is a short-range term, thus we have 
obtained an integrable bound h{t) for the integrand in (57). The integral exists as 
a Bochner integral or as an improper Riemann integral. Thus the limit exists for a 
dense set of states ip, and the wave operators exist as a strong limit on H. (We have 
not used the fact that diAk decays integrably, but it will be needed in the relativistic 
case [19].) In an arbitrary gauge A' for the given B, existence of the wave operators 
follows now from the transformation formula (65). The scattering operator is defined 
by S := Q*^Q_ . The Schrodinger equation is treated analogously, by omitting the 
term — cr-B. (Existence of Q± can be shown for every medium-range A, without any 
assumptions on div A or curl A, by quadratic form techniques. This proof requires 
an additional regularization, and it is not suitable for obtaining a high-energy limit.) 



4.3 Gauge Transformation 

Suppose A and A' are vector potentials of medium range with curl A = curl A' in 
S', thus A' = A + grad A and A(x) = limr_,oo A(rx) is continuous on W \ {0} by 
Thm. 1.2. We claim 

s-lime^^o^A(x) e^^^^t = A(±p) . (61) 

t— >±oo 
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Note that (^A(x) — A(x)j {Hq + i) ^ is compact, thus 

s-hmei^o^ (A(x) - A(x)) e^^^o^ = . (62) 
Moreover, since A is 0-homogeneous, we may multiply the argument with ±m/t > 0: 

e^^o^A(x) e~^^o^ = A(x + tp/m) = A(±mx/t ± p) (63) 
= e-i"^^V(2t)A(±p)ei^^V(2t) ^ A(±p) (64) 

strongly as t — + ±oo, since x^/t pointwise for x G M*^. This proves (61). Now 
consider the Hamiltonians H and H' = e^"^^^^ He~^'^^^^ according to Lemma 4.1, 
and suppose that Q± exist. Then the wave operators 

n'^ := s-hmei^'^e-i^ot = s-lime^^W (e^^^e-i^ot) (eii^ot,-iA(x)^-ii7ot) 

= ei^Wfi±e-i^(±P) (65) 

exist as well, and the scattering operators satisfy S' = e^'^^^^ S e~^'^^~^\ The 
gauge transformation formula is employed, e.g., in [31, 35, 37, 40]. The proof (64) 
seems to be new. The analogous formula for the Dirac equation is found in [17, 18, 
19]. If A — A' is short-range, then A is constant, and S' = S. 



4.4 Asymptotic Completeness of Wave Operators 

The wave operators Q± are called asymptotically complete, if every "scatter- 
ing state" in the continuous subspace of H is asymptotic to a free state, i.e., 
Ran(fi_) = Ran(fi+) = H™"*(i7) = W%H), which implies that S is unitary. Con- 
sider again the special gauge A = A* + A*": In the case of A'* = and Aq = 0, 
completeness was shown in [21, 7] by the Enss geometric method. The proof shall 
extend to our case, since the additional short-range terms can be included with 
standard techniques, but I have not checked the details. In an arbitrary gauge 
A' = A + grad A, completeness is carried over by the gauge transformation (65): we 
have Ran(fi'±) = e^^^^) Ran(f]±) = e^^'^^'>n'"'{H) = W{H'). Under the stronger 
assumptions B G Lf^^ and curlB(x) = C(|x|"*^^+'^^), completeness was shown by 
Arians [1] in the transversal gauge. Her proof employs a phase space cutoff of the 
form /(p — A(x)), which can be defined by a Fourier transform. 



4.5 Modified Wave Operators 

For A of medium decay, the unmodified wave operators exist although A(x) need 
not decay integrably. We shall compare them to modified wave operators: These 
exist when A(x) = A'^(x) + A'(x), where A'^ is short-range, and A' is C°°, with 
A'(x) = Odxl"*^), and with decay assumptions on its derivatives. The Dollard 
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wave operators are obtained from a time-dependent modification. If A is of medium 
range, it could be done in the form 

:= s-lime^^^[/^(t) , (t) := exp { - iHot + i f' p ■ A\sp) ds} , (66) 

thus = fi-t exp |i J^°° p ■ A'(sp) c/sj. By choosing A' transversal, the modifica- 
tion is vanishing. The time-independent modification of Isozaki-Kitada is employed, 
e.g., in [35, 25, 30, 39]. With Fourier integral operators J± we have 

ni := s-lime^^^ J± e''^^^ , J± : e^^^ ^ m^(x) (67) 

i— utoo 

for smooth A(x) = Odxl^"*). Here 'U±(x) is an approximate generalized eigenfunc- 
tion of H with incoming or outgoing momentum q, e.g., according to [39]: 

M±(x) := exp I iq • X + i y uj> ■ A{ujs) — uj ■ A(x + ujs) (is | , q = |q|uj . (68) 

Under a change of gauge, A — A' = A + gradA with grad A(x) = (9(|x|~^), the 
FIOs and the modified wave operators are transformed according to 

= e^-^ W - J± , l]f = ei^ W - , (69) 

and the modified scattering operator := '^'^j^ is gauge-invariant. (This is not 
the case when a; • A(a;s) is omitted from the integrand in (68), or in the idealized 
Aharanov-Bohm experiment, where A(x) is unbounded at x = and A(x) = A(x) 
[25, 37, 30]. Then Q'l^ = Q± and S'^ = S transform according to (6).) When A is 
smooth and medium-range, [24, Lemma 2.2] implies 

gj ^ g-ia(p) 5ei«(-P) , a(p) := Hp- A{sp) ds . (70) 

JO 

This may be taken as the definition of when A is not smooth. Cf. Sec. 6.2. 



5 High-Energy Limit and Inverse Scattering 

Consider the scattering of a state ^^^ip^ where u = uu. The position operator x is 
generating a translation by u in momentum space, and the limit of e 

-iux^^giux^^) 

gives the high-energy asymptotics of the scattering process, as m — oo for a fixed 
direction uj G S^~^. 



5.1 High-Energy Limit of S 

Applying the momentum-space translation by u = uuj to the free Hamiltonian gives 

e-i"^i/oei"^ = -^(p + v.f = ^u' + -L:.p + ^Ho) . (71) 
2m 2m m^ u ^ 

In the free time evolution, the first term is a rapidly oscillating phase factor, which 
cancels with the corresponding term for H in the Cook integral (72). Rescaling 
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the time t' = ut/m, the free time evolution is generated by a; ■ p + ^Hq ^ a; • p 
as n ^ oo. For the Hamiltonian H, note that A ■ p becomes A ■ (p + u) before 
rescaling, and H is replaced with if^ := a; ■ (p — A) + ^H. Employing the special 
gauge A = A* + A'' according to item 4 of Thm 1.2, consider 

e-i"^fi±ei"^^ (72) 

= ^ + i e^Hut {-u;.A + ^{H-Ho)}e-^(^P + ^H,)t ^ 
Jo ^ u ^ 

Assume ip & S with 'ip G . The velocity operator corresponding to the translated 
and rescaled time evolution is a; + p/w. Fix < e < 1 and Uq > 0, such that 
supp-^ is contained in the ball |p| < Mo(1 — e), then the speed is bounded below 
by e for u > Uq . By the standard techniques from Sec. 4.2, i.e., the decomposition 
(59)-(60), an integrable bound h{t) is obtained for the integrand in (72), uniformly 
for u > uq . The critical term is — A*" ■ (a; + p/w) = — G ■ [x x (a; + p/u)] with 
A^ = G X X. Again, the translated angular momentum x x (a; + p/n) is commuting 
with the translated free time evolution, and G(x) is short-range. By the Dominated 
Convergence Theorem (for the 7Y- valued Bochner integral), the limit m ^ oo is 
interchanged with the integration: 

lime-i"^l]±ei"^^ = ^ + i ei'^(P " (-a; ■ A) e'^^P^ rft (73) 



lim ei^(P- A)te-i^P^^ (74) 

r±cxD 



= exp|-i^ ■ A(x + a;t)rft}?/; . (75) 

We have employed the fact that ^ a; ■ (p — A) in the strong resolvent sense. 
The last step is verified from a differential equation [2, 18]. A density argument 
yields strong convergence of VL± , and (7) follows from S = and the strong 

convergence of Q'^ . (If Qu are isometric, Q^o is unitary, and Qu ^oo strongly, then 
Jl* — >• strongly because ^2* — fij^ = f2* (i^oo — ^u)^lo •) arbitrary medium- 

range gauge A' consider A' = A + grad A , the limit A according to Thm. 1.2, and 
the transformation formula (6). We obtain 

s-lime-i"^5'ei"^ = e^^H exp |i T a; ■ A(x + a;t) dtj e'^^^-^^) 

= exp |i /" (jj ■ A'(x + (jjt) dt\ , 

^ J — oo ^ 

since A(x + ujt) — > A(±a;) pointwise for t ±oo. Thus the high-energy limit (7) 
is established for an arbitrary gauge A. Under stronger decay assumptions, error 
bounds are given in Cor. 5.3. The same proof applies to the Schrodinger equation 
by omitting the term — cr ■ B. 

The simplification of the scattering process to a mere phase change at high energies 
has a geometric interpretation, cf. [8, 10]: For a state ijj with momentum support in 
a ball of radius mR around u = uu: and large u = |u|, translation dominates over 
spreading of the wave packet in the free time evolution. On the physical time scale 
t, the region of strong interaction is traveled in a time t ~ m/u, and the effective 
diameter of the wave packet is increased by ~ -Rt ~ 1/m. 
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5.2 Reconstruction of B 



In the inverse scattering problem, Aq , B and A are unknown, and the high-energy 
hmit of S is known up to a gauge transformation. The absolute phase of (7) is not 
gauge- invariant, but we assume that the relative phase is observable. The exponen- 
tial function is 27ri-periodic, but by the continuity of A, the integral transform of 
A is obtained up to a direction-dependent constant, and the magnetic field B is 
reconstructed according to Prop 3.1. This concludes the proof of Thm. 1.4 for the 
Pauli operator and the Schrodinger operator. 

5.3 Reconstruction of 

For the Schrodinger equation with short-range electromagnetic fields, Arians [2] first 
reconstructs B from (7), and then she reconstructs Aq from the l/u-term of the high- 
energy asymptotics. The following theorem is quite similar, but the proof will be 
modified. Similar results for B e C°° are obtained by Nicoleau in [24] using Fourier 
integral operators. 

Theorem 5.1 (Reconstruction of Aq (Arians)) 

Consider a short-range electrostatic potential Aq according to (4), and a magnetic 
field B : M'^ ^ M'^ or B : ^ M. Suppose B and curl B are continuous and both 
decay as 0{\x.\~^) with /i > 2, and that the flux o/B is vanishing if u = 2. A is a 
corresponding vector potential of short range. Set a(a;, x) := u: ■ A(x + ujt) dt. 
The scattering operator S for the Schrodinger- or Pauli equations has the following 
high-energy asymptotics: 

/oo 
Aq{^ + ujt) %l) dt 
-oo 

ie^^ r K^!'(x + a;t, p)7/^rft (76) 

J —oo 
roo 

i / K'^(^ + ujt, p)e^^i:dt 
Jo 

fori/j G . The operators K^{'x, p) are defined in (81), they depend on B but not 
on Aq. Thus B and Aq are reconstructed uniquely from S . 

Proof: For every a; G S'^~^ and the signs ±, consider the vector potentials 

/•±oo 

A^(x) := / xB{^ + ujs)ds = A(x) + grad A^(x) (77) 
Jo 

r±oo 

with A+(x) := / uj ■ A{x + u}s)ds . (78) 

Jo 

This definition is motivated by [35, 2, 24], and the transformation (77) is verified with 
a; X (V X A) = V(a; ■ A) — (a; ■ V) A. A^ and grad are bounded and continuous. 
A^(x) decays as |x|~*^^~-^'' in the half-space a; ■ x — ±oo, but in general it does not 
decay for a; ■ x ^ =Foo, thus A^ is not of medium range. We have curl A^ = B 
and div A^(x) = — J^°° u ■ curlB(x + a;s) ds in S'. The latter function is bounded. 



lini ^e-i"^f^-ei°)ei"^^ = - 

Li — >oo 777 \ / 
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continuous, and decays as |x| for cj-x ^ ±00. The Hamiltonian is defined 

by (54) with instead of A, it satisfies = e^'^i^^) H e'^^^^^^^ Analogously 
to the gauge transformation formula (65) we have 

f]^^ = e-i^±W lim e^^± ^ e~^^o^ , (79) 

t — ^±00 

when the support of is contained in the half-space ±0; ■ p > 0. To verify that 
the analog of A(±p) is vanishing, a nonstationary phase estimate is employed for a 
dense set of states. For ip G C^, consider the translated and rescaled Cook integral 

^g-iux (eiA^(x) _ ^\ giux^ (gQ) 
= i jT^^ ^K^P + ^H^)t {^^(x) + ir-(x, p)} e-H'^P + fHo)t ^ 
with the symmetric operators (omit —cr ■ B in the Schrodinger case) 

i^^(x, p) := ^ ( - 2A^(x) ■ p + idiv A^(x) + (A^(x))2 - a ■ B(x)) . (81) 

Note that uj ■ A^ = 0, thus the main contribution from (72) is missing, and a 
common factor 1/u was extracted from the remaining terms. By the same uniform 
estimates as in Sec. 5.1, the limit n ^ 00 can be performed under the integral. Since 
the translated and rescaled generators of the time evolutions are converging to ■ p 
in the strong resolvent sense, (80) is converging to 

[Ao{x + ujt) + K^ix + Ljt, p)}ijdt . (82) 

The relation a(a;, x) = A':f'(x) — A'^(x) gives (76) as a weak limit. Moreover, we have 

e"^^''^ e^^"^ — >■ e^'^i^''^-' strongly, and the strong limit in (76) is obtained from 
the decomposition 

S -e^^ = Ql - e-i^'^) + fi; (e~i^+ - Q+) e^^ . (83) 

For the limit of f2+ , we are applying (80) to e^*^ ip G W'^ instead of ip, which re- 
quires two additional arguments: First, the regularization of Aq is not done with 
^ = g{p)^lj, but with e^" = {Ho + i)~^ ({Hq + i) e^^ i?j . Second, in the proof of the 

nonpropagation property, {e^'^ ip)^ does not have the desired compact support. But 
a{iv, x) is constant in the direction w, thus the support of ip is enlarged only or- 
thogonal to uj, and remains bounded in the direction of cj. (Use only the directional 
derivative u: ■ Vp in the proof of the nonstationary phase estimate [36, Thm. 1.8].) 

Now suppose that 5* is known (it is invariant under short-range gauge transfor- 
mations), thus the absolute phase of its high-energy limit is known. Then B is 
reconstructed by Thm. 1.4, and any corresponding short-range A gives a{u}, x). 
Since -R'±(x, p) can be computed from B, the X-ray transform of Aq is obtained 
from (76). Aq is reconstructed in the second step according to [14], at least under 
stronger regularity assumptions. In the general case, the potential is regularized by 
translated test functions [9], or it is considered in S' [18]. H 
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The main difference to Arians' original proof [2] is the adaptive gauge transformation 
A" = A + grad A" : since u ■ A^(x) = 0, the limit (82) is read off easily from (80) 
after showing the uniform bound, and the expression (81) for A'^(x, p) is obtained 
from if" — Hq without calculation. In [3], magnetic fields with compact support 
and nonvanishing flux are considered as well, by employing a family of transversal 
gauges with adapted reference point. Cf. Cor. 2.8. Analogously we have 

Remark 5.2 (Generalization (Arians)) 

Suppose V = 2 and , B satisfy the assumptions of Thm. 5.1, except the flux is 
not vanishing. Then (76) and the proof remain valid, if A(x) decays integrably in 
the half-plane a; ■ x > and in a sector around —uj. When A is fixed, we may 
consider a family of gauge transformations according to Cor. 2.8 to satisfy the decay 
requirements, and the right hand side of (76) is modified. 

As noted in [2, 3], the uniform estimate of the Cook integral (80) and Remark 5.2 
give error bounds for the high-energy limit of S according to Thm. 1.4: 

Corollary 5.3 (Error Bounds (Arians)) 

1. Under the short-range assumptions of Thm. 5.1, the limit (7) has an explicit error 
hound for G , which is of the form 

^-iux^^iux^ = ei«(^' ^) ^ + 0{l/u) . (84) 

2. If u = 2 and the flux of B is not vanishing, then an analogous weak estimate 
remains valid for <f), ip E and medium-range A with A(x) = (9(l/|x|).- 

(0, e-i"^^ei"^V^) = (0, ei«(^' ^) ^) + 0(l/w) . (85) 

Error bounds for (76) would require stronger decay assumptions on Aq [9], or 
stronger regularity assumptions [24]. The right hand side of (76) contains a multi- 
plication operator times p. In the short-range case, it can be rewritten according 
to 

/O poo POO 

A'^(x+a;t) dt+ / A^;^(x+a;t) dt = ta;xB(x+a;t) dt = a{cj, x) . (86) 
~oo Jo J —oo 

Only the last expression remains meaningful, if B is smooth but does not not decay 
faster than |x|~^: an asymptotic expansion of e~''^'^ S'e^^'^ in powers of \/u has 
been obtained in [24] for B e a"B(x) = C(|x|-^-l"l), /i > 3/2. 

In the special case of B = and A = 0, thus H = Hq + Aq the uniform estimate 
of (80) together with S - I = - fi+) and the strong limit of Vt\ yield 

lim - f e"^"^ S e^"^ -l)^ = -i Ao(x + ujt) ^ dt (87) 

u^oo m ^ ' J-oo 



for ip G C^. This limit was obtained by Enss and Weder in a series of papers, cf. [8, 
9], introducing the time- dependent geometric method, and including also the cases 
of A^-particle scattering, inverse two cluster scattering, and long-range electrostatic 
potentials. Error bounds for the weak formulation of (87) were established in [9] 
under stronger decay assumptions. The strong limit is due to [17]. 
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6 Concluding Remarks 



A geometric interpretation of the scattering process at high energies is given in 
[8, 10], cf. Sec. 5.1. Here we shall discuss the implications of gauge invariance and 
the question of measurable quantities in a scattering process, as well as the possible 
application of the inverse scattering problem. 

6.1 Gauge Invariance 

The vector potential A is determined only up to a gradient by the magnetic field 
B (and by additional discrete values of fluxes, when the domain is multiply con- 
nected). It is not eliminated easily from the Schrodinger equation iip = Hip. (For 
the nonlinear hydrodynamic formalism, cf. [33, 31, 26].) The Schrodinger equation 
or Pauli equation is invariant under the simultaneous gauge transformation of A 
and ip, 

A ^ A' = A + grad A tp ^ ip' = e^^ . (88) 

(If the electromagnetic field was time- dependent, we would have A'q = Ao — X in 
addition.) One interpretation is, that the transformation of A comes from curl A = 
B, and it needs to be compensated for by the transformation of ip. Following Weyl, 
this argument can be reversed: if we assume that the local phase of the wave function 
ip is not observable, the theory should be invariant under a local U{1) transformation, 
and Ip must be coupled to a vector potential. It is assumed in general that A and 
A' describe the same magnetic field, and that all observable physical effects are 
independent of the chosen gauge. Thus an electron in an electromagnetic field is 
described by an equivalence class of pairs (A, tp), where (A, tp) ~ (A', tp') iff there 
is a A with (88). 

Gauge invariance implies that not every self-adjoint operator corresponds to a phys- 
ical observable: Suppose that the self-adjoint operator F = F(p, A(x), . . .) is con- 
structed from a function /(p, A(x), . . .) by some quantization procedure (since the 
ordinary functional calculus does not apply due to [x, p] 7^ 0). Then the expecta- 
tion value {ip^ Fip) is gauge- invariant, i.e., it depends only on the equivalence class 
of (A, ^), iff 

gi^W F{p, A(x), . . .) e~i^W = F(p, A(x) + grad A(x), . . .) . (89) 

This restriction is a "superselection rule" in the general sense of [33]. At least if / 
is polynomial in p, this means that it depends on p and A only in the combination 
p — A. Examples of nonobservable operators are A(x), the canonical momentum 
operator p, the free Hamiltonian Hq = 2^p^, and the canonical angular momentum 
L = X X p. The following operators are among the observables: x, Aq[x.), B(x), the 
kinetic momentum mx = p — A, the kinetic energy — A)^, the Hamiltonian 

H, and the kinetic angular momentum x x mx = xx(p — A) = L — xxA. See 
[33, 19] for a discussion of vector potentials and gauge invariance in the context of 
the Aharanov-Bohm effect [26]. 
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6.2 The Scattering Cross Section 



When B = 0, or for the time evolution of asymptotic configurations, it is natural 
to set A = by convention. In the scattering theory with medium-range vector 
potentials, we must assume that A and A' describe the same physical system, as 
soon as curl A = curl A'. Then A' — A need not be short-range, and the scattering 
operator S is not gauge- invariant, but it transforms according to (6). Given a 
scattering state ^ Ran(f2+), the particle is found in a cone C for t +oo with 
probability 



lim 

t— >+oo 



F(x e C) e-i^* ij 



(90) 



according to Bollard [28, Thm. IX.31]. Now VL%iIj' = e^^(P) fi^?/; by (6), thus 
this number is gauge-invariant. By the correspondence between subsets of S'^~^ 
and cones in (with apex 0), (90) defines a measure on S"^~^. The differential 
cross section da/duj for incident momentum q is obtained when (p = is ap- 

proaching a plane wave, rescaled such that its Fourier transform 0(p) is approach- 
ing "•y/(5(pi — q) S{p2) ^(ps)" when q = (g, 0, O)*'". If the momentum support of an 
incoming asymptotic configuration is concentrated at p ~ q, or at p G q[0, oo), 
we have 

S'cp = ei^(P) 5e-i^(-P) ^ e'^(P) ^e-^^^-^) = e^^^P) " i^^'^) 50 . (91) 

The phase factor in momentum space does not influence the probability of finding the 
outgoing state in a cone, thus da/du is gauge-invariant. By the same argument, i.e., 
replacing A with a, we may compute it from the gauge-invariant modified scattering 
operator , which was defined in (70). See also [35, 30]. 



6.3 The Phase of the Scattering AmpUtude 

In short-range scattering theory, the differential cross section is — = /q(w) 
where fqiuj) is the scattering amplitude for incident momentum q and outgoing 
direction uj. It is obtained from the T-matrix, i.e., the integral kernel of S* — 1 on 
the energy shell. Probably this relation remains valid in the medium-range situation, 
although the latter kernel will be more singular on the diagonal [39]. In scattering 
experiments, usually the differential cross section is observed for an incident beam 
of particles, which is modeled as a plane wave. Information on the phase of the 
scattering amplitude is not available directly, but it is required for solving the inverse 
scattering problem, e.g., with (87). 

2 

For a small, central, scalar potential in M^, /q can be reconstructed from /q(<^) 
by employing the unitarity of S, see [27, Sec. V.6.B] and the references in [20]. If 
this approach is extended to R^, it will work for rotationally symmetric B as well. 
Phase information would be available experimentally, if it was possible to localize the 
incoming particles more precisely. It could be reconstructed as well, if the location of 
the unknown scatterer is kept fixed, and the location of a known additional potential 
is varied while measuring the cross sections [20]. In some cases, phase information 
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is obtained from interference between the scattered beam and a coherent reference 
beam [26]. 

6.4 Two-Particle Scattering 

If two nonrelativistic particles are interacting via a pair potential Ao(x2 — xi), their 
relative motion is equivalent to an external field problem for one particle with the 
reduced mass m = "^^.'"^ . The pair potential has a physical iustification only if it 
is a central potential, or if, say, m2 -C mi : Suppose particle 1 is a molecule with a 
dipole field given by Ao(x2 — xi), and particle 2 is an electron. The molecule will 
be rotated by interacting with the electron, but the corresponding rotation of Aq is 
neglected in the model. This simplification is justified if mi ^ m2 , and in this case 
we might assume as well that the molecule is generating a magnetic field. If the 
orientation of particle 1 is unknown, the scattering cross section may be defined by 
averaging over the orientations. But the phase information required for solving the 
inverse problem is unlikely to be recovered. 

6.5 Inverse Scattering and Error Bounds 

The high-energy asymptotics (7), (76), (87) might be applied independently from 
inverse scattering. But consider the inverse scattering problem for a particle in 
an external electromagnetic field, or for two-particle scattering under the restric- 
tions from Sec. 6.4. The uniqueness of the solution is interesting from a theoretical 
point of view, because in the analogous situation of particle physics, the models are 
based mainly on scattering data. In our situation, two additional questions must 
be addressed, before a field can be reconstructed from a scattering experiment: the 
problem of obtaining phase information, cf. Sec. 6.3, and the effects of a high but 
finite energy. When a-priori bounds on , B, and A are given, it is possible to 
estimate the approximation error in (7) according to Cor. 5.3. Thus we can check 
in principle, if the required energy is available in the experimental setup, and if the 
nonrelativistic model makes sense for this high energy (note also that the scattering 
operators for the Pauli- and Dirac equations (positive energy) coincide if = 0, 
cf. [36, 19]). 

If the high-energy limit was observed for a suitable family of states ip, we could 
obtain the required X-ray transforms as multiplication operators. At a high but 
finite energy, we do not get the operator of multiplication with an approximate X- 
ray transform, and it is not clear how to obtain the latter. (If the X-ray transform 
was obtained approximately, we could apply regularization techniques to provide an 
approximate inversion of the X-ray transform, whose exact inversion is ill-posed.) 
For Aq of compact support, it is suggested in [34] to consider (87) for a single chosen 
ip. Or we might specify an inversion procedure for the X-ray transform and apply it 
to the high-energy asymptotics, to check if the resulting composition of operators is 
converging. It may be possible as well, to obtain Aq at a lower energy by a recursive 
approach, or by considering more terms of the asymptotic expansion. 
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Corresponding Results on Relativistic Scattering 

This appendix of the preprint will not be part of the published paper. It should be 
considered as a summarizing preview of [19]. 

The Dirac operators are matrix-valued operators of the form Hq = a ■ p + (3m and 
H = a - {p — A) + (3m + Aq . The vector potential A shall be of medium range. The 
scalar potential shall be continuous except for a finite number of local singularities, 
where |y4o(x)| < c|x — Xj|~^ with fi < 1, and decay integrably. The scattering 
operator S is decomposed according to the subspaces of positive or negative energy, 
and we will need the Newton- Wigner position operator x [36]: 

Theorem Al (High-Energy Asymptotics and Inverse Scattering) 

Suppose that B a magnetic field of medium decay in or M.^, A is any medium- 
range vector potential with curl A = B m S' , and Aq is a short-range electrostatic 
potential. 

1. The wave operators fl± and the scattering operator S for the Dirac equation exist. 
Consider also a gauge transformation A' = A + gradA and A(x) = lim^^oo •^('"x), 
and denote the operators corresponding to A' by H' , Q'^, S' . They obey the gauge 
transformation formula 

n'^ = ei^W f]±e-i^(±P sign(^ro)) s'^ = ei^(±P) S±e-'^(^P) , (Al) 
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where S± denotes the restriction of S to the subspace of positive/negative energy. 

2. Consider translations in momentum space by u = uu), uj G S^~^ . Denote by 
S± the restriction of S to the subspace P±H of positive/negative energy. With the 
Newton-Wigner position operator St we have the high-energy limit 

s-lim e~^"^ S± e^"^ = exp { - i T (Aq t ■ a) (5 ± ut) dt} . (A2) 

Both Aq and B are reconstructed from the relative phase of the high-energy limit of 
S+. Analogous results hold for the Klein-Gordon equation. 

Sketch of the proof: 1. The Cook integral is estimated by employing the special 
gauge A = A'' + A'' and applying the technique of [22] to A^. We have A^ ■ a = 
■ p/Hq -\- ■ {ex. — p/Ho). The first term is written as A** ■ p = G ■ L, and 
the second term is controlled with partial integration, since it is oscillating due to 
{ex — p / Hq , Hq} = 0. The usual density argument gives the existence of fl± . The 
gauge transformation formula was proved in [17]. 

2. The high-energy limit of [18] is extended in two directions: Medium-range vector 
potentials are included by estimating the time evolution analogously to the direct 
problem. Local singularities of Aq are treated with a density argument and variable 
cutoff functions. 

The corresponding results for the Klein-Gordon equation are obtained analo- 
gously by writing it as a first-order system in the Foldy-Wouthuysen represen- 
tation and by employing the Dyson expansion. Now Aq shall be bounded and 
\\Ao{x)\/p^^rm^ II < 1. The proof was given in [10] under more restrictive condi- 
tions, it is extended again with A^ ■ p = G ■ L. ■ 

Consider a compact set K C M*^ and a Hamiltonian Hk in //^(M"^ \ K). Obstacle 
scattering means to deal with the wave operators 

n± := s-lime^^^^Je-'^o^ , (A3) 

t— >±oo 

where J is the natural projection from L'^{W^) onto L^(M'^ \ K). Nonrelativistic 
inverse obstacle scattering using high-energy limits was discussed by Nicoleau [25] 
and Weder [37]. We shall consider the Dirac equation in two cases: 

i) The compact set K C M'^ is convex and space-reflection symmetric, i.e., K = —K. 
a-p is a symmetric operator on C^{W\K), and the deflciency indices of its closure 
are equal, since it is anticommuting with the unitary involution {R%l)){x) = ?/;(— x). 
Fix any self-adjoint extension and employ the Kato-Rellich Theorem to deflne '■ = 
a ■ p + (/5m — a. ■ A + Aq) when Aq is bounded. 

ii) In the case of i^' = {0} C M^, stronger singularities of A at x = are permitted: 
the magnetic fleld shall be of the form B(x) = Bm(x) + $*5(x) with B^ of medium 
decay and G M. The vector potential has a decomposition A(x) = Am(x) + 
A*(x), where A^ is of medium range and curl A^ = B^. 

A.(x):=|:^|x|-f-"^) (A4) 
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satisfies curl A* = $^5, and it behaves as a vector potential of medium range for 
|x| oo. Now A^,(x) is odd, thus a ■ [p — A*) is a symmetric operator with equal 
deficiency indices. Fix any self-adjoint extension and include j3m — ct ■ A,„ + Aq with 
Kato-Rellich. 

Theorem A2 (Obstacle Scattering) 

Consider an obstacle K C W and a Dirac operator Hk satisfying Assumption i) 
or ii) above. 

1. The wave operators (A3) exist and are isometric, they transform under a change 
of gauge according to (Al). 

2. For uj G S^~^ and all V' G L^(R'^, C^) with ^-support outside of the cylinder 
K + M.(jj, i.e., -F(x G A' + M.uj)il) = 0, we have the high-energy limit 

w-hm e~^"^ S+ e^"^ i; = exp { - i T (Ao - ■ A){^ + ujt) dt] i' . (A5) 

3. For^ & My\K, the electrostatic potential AoijK.) and the magnetic field Blx.) are 
reconstructed from the relative phase in (A5). If u = 2 and K ^ {0}, we need the 
additional assumption that both decay faster than any power. 

4. If u = 2 and B has a finite flux $, then $ is reconstructed modulo 2n from the 
relative phase in (A5). 

If the absolute phase was observable, then $ could be reconstructed modulo 47r 
[25, 37]. Note that the obstacle is assumed to be known, and only the fields are 
reconstructed. Different self-adjoint choices of Hx are not distinguished in the 
high-energy limit. The stronger decay assumptions in item 3 are required by the 
Support Theorem for the X-ray transform [14]. In the case ii), items 3 and 4 mean 
that Bm is reconstructed uniquely, and is reconstructed modulo 27r. If z/ = 2 
and the flux of B on K is nonzero, it is influencing the particle outside of K via the 
vector potential (Aharanov-Bohm effect). 

Sketch of the proof: The Cook integral is estimated analogously to the case 
without obstacle, by replacing the projector J with a smooth cutoff function 
When interchanging the limits m — > oo and t — > ±c)0, the free time evolution is 
estimated by introducing a momentum cutoff f{p/u) as in [37], since the states do 
not have compact momentum support. For the high-energy limit at a finite time, the 
Dyson expansion does not apply when K ^ {0}, since J is not injective. Introducing 
the Dirac operator H in L'^{M.'^), consider the decomposition 

^iHKt^^-iHot^ = (ei^^^xe"^^^) (e^^^e~i^o^) (A6) 

= X e'^^ e-i^o^ ^ + i /* e^^^^ (HkX - xH) e'^^^ ' ^) e-^^o^ ds . (A7) 

Jo 

The high-energy asymptotics of the time evolution are known for Hq and H but not 
for Hk. After performing the known limits, the integral is seen to vanish because of 
the support properties of ip and and the limit of the first term in (A7) remains. ■ 
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